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I. INTRODUCTION 



Avoided crossings 0| have been observed theoretically and experimentally in a large 
variety ofphysical systems involving eigerivalues (e.g., natural frequencies, eigenenergies) 
0, S B B B 0, i, i, H Q Q 13 13, Q, and they have attracted much attention 
even in recent years because of their rich physics and practical importance in, for example, 
mechanical engineering [l^, 13 H| and quantum physics 0, |E 0, B ■ avoided 
crossing regions, eigenvalues of the system first approach each other as a system parameter 
is varied but then veer abruptly from each other without crossing. In those regions a drastic 
change of some characteristic of the system occurs along the eigenvalue loci. In Ref. j3, fo^^ 
example, Pierre illustrated that the mode shapes of a disordered chain of coupled pendulums 
change in the regions where avoided crossings of the eigenfrequencies of the system take 
place. In that study, disorders in the lengths of the pendulums were used as the system 
parameters. Also, in Ref. 0, Walkup et al. studied in detail avoided crossings observed in 
the energy levels of diamagnetic hydrogen as functions of the magnetic field strength or the 
angular momentum, which lead to the diabatic exchange of the states of the wave functions. 
A study of Bose-Einstein condensation (BEC) B showed that in order to trap molecules 
created in an atomic BEC through a Feshbach resonance, an avoided crossing of two bound 
states of the molecules must be exploited, through which the vibrational quantum number 
and size of the trapped molecules change. 

In the present paper, we show theoretically that avoided crossings can be observed in 
acoustically coupled bubbles, which has to the authors' knowledge not been stated in the 
literature. Furthermore, based on analyses of transition frequencies we propose a 

way to detect a state exchange occurring in the avoided crossing region. The theoretical 
model used in this study, reviewed in Sec. HIl is a forced coupled oscillator model that 
describes acoustic coupling of pulsating bubbles. Using the model, we show in Sec. Illll that 
if at least three bubbles exist, it is possible that the resonance frequencies of the bubbles 
exhibit an avoided crossing when they are plotted as functions of the separation distances 
between the bubbles. As has been demonstrated (e.g., Refs. [3, 13); double-bubble 
systems, neither crossings nor avoided crossings of the resonance frequencies as functions 
of the separation distance are observed, since the higher of the two resonance frequencies 
of the systems increases and the lower one decreases as the bubbles approach each other. 
However, as shown in the present paper, by introducing one more bubble whose monopole 
(i.e., decoupled) resonance frequency crosses with one of the resonance frequencies of a 
double-bubble system, one can observe the avoided crossing of the resonance frequencies 
when all three bubbles are coupled. 

In Sec. IIVI we examine the phase properties of the three coupled bubbles to show that 
a state exchange actually occurred between the bubbles in the avoided crossing region. In 
this effort, the notion of a transition frequency plays an important role. The transition 



frequencies introduced in Refs. [16|, ll7| are characteristic frequencies of acoustically coupled 



bubbles, around which the oscillation phase of bubbles inverts, e.g., from in-phase to out-of- 
phase with the driving sound. It was proved in Ref. ^3] that a bubble in a A^-bubble system 
has up to 2N — 1 transition frequencies, only N ones of which correspond to the resonance 
frequencies of the system. That is, observing the transition frequencies allows us to obtain 
richer insight into the phase properties than that obtained by only observing the resonance 
frequencies. This notion has already been exploited as a powerful tool to understand the 



-Pioitea 

sign reversal of the secondary Bjerknes force 0, in which the oscillation phases play a 
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crucial role. Using this notion and observing directly the oscillation phases, we show that 
the coupled bubbles exchange their oscillation states through the avoided crossing and the 
state exchange takes place at the separation distances where an avoided crossing resonance 
frequency crosses with a transition frequency that is not a resonance frequency. The present 
findings appear to reveal a taste of bubbles' hidden complexity. 

Section |3 summarizes this paper, and the Appendixes present additional remarks. 



II. COUPLED OSCILLATOR MODEL, RESONANCE FREQUENCY, AND 
TRANSITION FREQUENCY 

The theoretical model used in the present study is a forced oscillator model in which 
harmonic oscillators are coupled and references therein): 

r, 1 ^ 7?2 

Ci + Ui^ei + diei = -—— - — ^ j^ej (1) 

prCio -KjO j^i J^ij 

for i = 1,2,---,N, 

where N corresponds to the number of bubbles, Rio is the equilibrium radius of bubble i, 
Ci is the deviation of radius assumed as |ej| <C Rio, Uio is the monopole (angular) resonance 
frequency of bubble i, defined as 

(2) 



\ pRl 



Si is the damping factor, the overdots denote the time derivation, is the pressure of 
the external sound, p is the density of the surrounding liquid, Dij (= Dji) is the separation 
distance between the centers of bubbles i and j, Hi is the polytropic exponent of the gas inside 
the bubbles, Pq is the static pressure, and a is the surface tension. In this linear model, 
the following assumptions are made: the surrounding liquid is incompressible, the sound 
amplitude is sufficiently low, the separation distances are much larger than the bubbles' 
radii, and the shape deformation of the bubbles is negligible. The last term of Eq. (P), 
representing the pressures of the sounds that the neighboring bubbles emit, describes the 
acoustic coupling between the bubbles. As in the double-bubble case jl^, this model may 
be assumed to be of third order with respect to the inverse of the separation distances (i.e., 
the truncated terms are of fourth or higher order); see Appendix 1X1 

Using this model with = 3, a matrix equation for determining the amplitudes and 
phases of the radial oscillations is derived. Assuming = —Pa exp(ia;t) and Cj = /3j exp(iu;t) 
with Pa being a positive constant, uj being the driving (angular) frequency, and /3j being a 
complex amplitude, we have 



P. 



/?2 = (3) 
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P 



where A is a 3 x 3 matrix whose elements, a^j- (i, j = 1, 2, 3), are defined as 

~ ' — —X otherwise, ^^"^ 

A, 
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with 

X = u;^ (5) 
and I = (1, 1, 1)"^. We should note here that essentially the same matrix equations can be 




found in previous papers (e.g. [17|, |2l|, |22| ) . The solution of Eq. Q is represented as 

= -jA-H (6) 

Pa \A\*CI 



P |A|*|A| 

where |A| and C are the determinant and the cofactor matrix of A, respectively, and |A|* 
is the complex conjugate of |A|. We used here an expression in which the denominator is 
real. 

The eigenfrequencies of the system are determined by 

|A| = 0, (7) 

which is a cubic equation in terms of X. For (5j ?a 0, the roots of this equation are equivalent 
to the resonance frequencies of the system. The transition frequencies of bubble i, defined 
as the driving frequencies at which the phase difference between bubble i and the driving 
sound is n/2 (or 37r/2) [3, 12) 13; determined by 



Mr^) = 0, (8) 

where 

T2\=\A\-Cl. (9) 

(See Appendix iBl for the concrete forms of | A| and C I.) From the mathematical proof given 
m Ref. hj, one knows that Eq. (jHJ is a fifth-order polynomial in terms of X, meaning that 
the bubbles may have up to five transition frequencies. 

The phase delay of bubble i, denoted by (pi, measured from the phase of the driving sound 
is determined using the atan2(?/,x) function in the C language, which returns tan~^{y/x) G 
[—71, 7i], as 

(k. = [^i if i)i > 0, 
\il>i + 2iT otherwise 

with 

ipi = atan2(— Im(rj), Re(ri)). 

The next section shows that in certain cases an avoided crossing is observed in the solution of 
Eq. ((7j). In the discussion, to obtain real eigenfrequencies that correspond to the resonance 
frequencies of the triple-bubble system for weak damping, we for the moment assume 5^ ~ 
(but 6i 7^ 0). Under this assumption, one obtains 

Im(|A|) ^ 0, (10) 

|A| ^ |A|*, (11) 

and 

r.^Re(r,). (12) 
Influences of the damping effect on the phase properties are briefly discussed in Sec. |TVl 
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III. AVOIDED CROSSINGS OF RESONANCE FREQUENCIES 



To begin with, a double-bubble system is briefly reconsidered to confirm that no avoided 
crossings are observed in the resonance frequencies of the system as functions of the separa- 
tion distance. The solid lines in Fig.nindicate the resonance frequencies of two coupled bub- 
bles (bubbles 1 and 2) of (-Rio, -R20) = (50 fim, 51 /xm) as functions of /12 = Di2/{Rio + -^20)- 
The other parameters are set to p = 1000 kg/m^, = 1.4 (i = 1, 2, 3), Pq = 1 atm, 
and a = 0.0728 N/m. As has been proved theoretically [ifl |20|, |23|, |2J , two resonance (or 
natural) frequencies appear in this system, each of which, for D12 00, converges to the 
monopole resonance frequency of a bubble. The higher resonance frequency increases and 
the lower decreases as the separation distance decreases. It is therefore obvious that avoided 
crossings cannot occur. 

Here we introduce one more bubble into the system. The dashed line displayed in Fig. ^ 
denotes the monopole resonance frequency of the introduced bubble, bubble 3, whose radius 
R30 = 51.5 fim. Note that this resonance frequency crosses with a resonance frequency of 
the double-bubble system. This crossing, as shown immediately, triggers an avoided crossing 
when the third bubble is coupled with the double-bubble system. 

Figure 121 shows the resonance frequencies in the case where all three bubbles are coupled. 
The separation distances are set to Du = lu^Rw + -^20), -D23 = h3{R2Q + R30), and D31 = 
D12 + D23] that is, the bubbles are arranged in line (see Fig. Ol^a)). Here the nondimensional 
quantities I12 and I23 are used as the system parameters. Figures Efa), |2tb), andEfc) show 
the results for Z23 = 100 [i^, 50, and 20, respectively. In the figures, an avoided crossing 
is clearly seen that takes place around the point at which the two decoupled resonance 
frequencies cross. The line of the resonance frequency originating with bubble 3 is divided 
into two parts, and each of them connects smoothly, like blending, with the curve of a 
resonance frequency of the double-bubble system, also divided into two parts. As bubble 
3 comes closer to the others, the avoided crossing becomes broader and the origin of each 
resonance frequency becomes increasingly unclear. 

An avoided crossing is also observed when bubble 3 is smaller than the others. Figure 
13] shows the resonance frequencies when R^q = 49.5 pm. Here the bubbles are arranged 
as illustrated in Fig. El^b). If bubble 3 is so large or so small that its monopole resonance 
frequency does not cross with a resonance frequency of the double-bubble system, no distinct 
avoided crossing is observed, though this situation is not shown here. 



IV. STATE EXCHANGE IN THE AVOIDED CROSSING REGION 

To manifest a state exchange like that which the coupled bubbles experience through the 
avoided crossing, we examined the oscillation phases of the bubbles. In bubble dynamics, 
the phase of radial oscillation plays important roles in many situations, including acoustic 
levitation ji^ 27, 2^, bubble-bubble interaction and multibubble sonoluminescence 



[29[, and hence an accurate understanding of it is crucial. In fact, by carefully examining 
the oscillation phases of two coupled bubbles for weak driving, we have recently succeeded 
in presenting a novel interpretation, which may be more accurate than previous ones, of 
the sign reversal of the secondary Bjerknes force 3 3; ^ paradoxical phenomenon that is 



considered to be the cause of the stable structure formation of bubbles in a weak acoustic 
field ,3fli] . In that discussion, it was suggested that the transition frequencies seem to be 
essential components for gaining an accurate understanding of the phenomenon, since the 
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sign reversal takes place at the transition frequencies that cannot be obtained by resonance- 
frequency analysis. In the present paper, we show by examining the oscillation phases 
that the bubbles exchange their oscillation states through the avoided crossing. As shown 
later, the point at which the state exchange occurs can be clearly detected by observing the 
transition frequencies. 

Figure shows the transition frequencies for (i?io, -R20; -R30) = (50 /im, 51 /zm, 51.5 /im) 
with I23 = 20. The thick lines denote the transition frequencies that correspond to the 
resonance frequencies already shown in Fig. |2fc). As expected from the mathematical proof 
presented in |l7[, the bubbles have up to five transition frequencies, all of which invert 
the oscillation phase of the corresponding bubble. It is worth noting that in each panel of 
Fig. El the second-highest resonance frequency (denoted below by ci;2nd) crosses once with 
a transition frequency in the avoided crossing region. Such crossings have not been found 
in double-bubble systems [l8| . In the following discussion, we focus our attention on 
the phase properties of the bubbles in this region to elucidate what happens around the 
intersecting points. 

The phase delays (pi for different lu as functions of u/ Uiq are shown by the solid lines in 
Fig. ini In the computation of (pi, we used very small but nonzero Si to obtain continuous 
results. Figures a, b) and [Hfc, d), respectively, show (pi for I12 smaller and larger than 
the intersecting point lu = hnt (~8.89). Here, we only displayed (pi in the frequency range 
around the two avoided crossing resonance frequencies. The vertical dotted lines indicate 
the two lowest resonance frequencies. As in double-bubble cases jl^ 3, at the resonance 
frequencies the phase delays of all bubbles shift simultaneously by whereas at the 
remaining transition frequencies only one phase delay shifts by — tt. 

The 0j-curves, as can be clearly seen in the figures, have different convexities on different 
sides of the intersecting point. For I12 smaller than Zmt, at a;2nd, 0i and 02 shift from vr to 
27r but (p3 shifts from to vr as increases. For /12 larger than Zjnt, on the other hand, an 
opposite tendency is seen; (pi and (p2 shift from to tt but (p^ shifts from tt to 27r. That 
is, a kind of state exchange takes place between bubble 3 and the other two bubbles at the 
intersecting point. 

Regarding the relationship between the state exchange and the phase properties, in the 
frequency range around ti;2nd5 bubble 3 oscillates out-of-phase with the other bubbles re- 
gardless of whether I12 < hnt or I12 > hnt, although the individual phase delays experience 
rapid shifts at ci;2nd- This means that the state exchange cannot be perceived accurately 
by observing whether the bubbles oscillate in-phase or out-of-phase with each other or by 
observing the sign of the secondary Bjerknes force, which is determined by the cosine of the 
phase difference between two bubbles ]31, i^l- Just the individual phase delays (or transition 
frequencies) should be examined. 

In the 0j-curves, we can find several similarities with double-bubble cases. Bubbles 1 and 
2, or bubble 3, have a phase delay greater than vr in the frequency range from ti;2nd to a 
certain higher frequency (equal to the next-higher transition frequency of the correspon ding 
bubble). A similar observation can be found for double-bubble systems 3 3- In Ref. |l8l | 
we discovered and elucidated that such a large phase delay can appear when two bubbles 
interact with each other through sound. In the double-bubble case, the larger one of the two 
bubbles has a phase delay greater than vr in the frequency range between the higher of two 
resonance frequencies and the highest of the transition frequencies of the bubble. We can, 
for a wider frequency range, also find a similarity between the double- and triple-bubble 
cases. In the frequency range uj/ujiq < 0.995, the profiles of (pi and (p2 for I12 < hnt and that 
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of 03 for li2 > lint are very similar to the profile of the phase delay of the larger bubble in 
a double-bubble system; those phase delays first exhibit two sharp rises and then one sharp 
fall as u increases. Also, the profiles of the remaining phase delays are very similar to that 
of the phase delay of the smaller bubble in a double-bubble system, exhibiting one sharp 
rise, one sharp fall, and then one sharp rise. This seems to indicate that in the frequency 
range considered, for /12 < hnt bubbles 1 and 2 act as "larger bubbles" while bubble 3 acts 
as a "smaller bubble," but for I12 > knt each bubble acts in the opposite way; that is, 
the physical roles that the bubbles play are exchanged through the avoided crossing. This 
observation could also be interpreted as a result of the change of a physical meaning of co'2nd- 
As illustrated in Fig. |21 ci;2nd is a hybrid of two resonance frequencies having different origins. 
We assume here that the origin, or the principal origin, of each avoided crossing resonance 
frequency is switched at Zint. This assumption allows us to consider that co'2nd for lu < knt, 
for example, is the resonance frequency whose principal origin is bubble 3. This suggestion 
is consistent, not only with the observation for large I23 where the origin of each resonance 
frequency is relatively clear, but also with the above speculation that bubbles 3 acts as a 
"smaller bubble" for I12 < knt, because a;2nd is higher than the lowest resonance frequency 
that is one of the two avoided crossing resonance frequencies. The observation for lu > knt 
can be interpreted in a similar manner. 

Lastly, we briefly examine how the damping affects the state exchange. For the damping 
coefficient, we use the value for viscous damping. 



4/. 



pR 



(13) 



with viscosity fi = 1.002 x 10~^ kg/ (m s). The dashed curves in Fig. IHlshow the phase delays 
in the damped case. The viscous effect smoothes the phase profiles, but the convexity of 
the curves is not altered from that for 5i ~ 0, as in the double-bubble cases [iiQ. The 
state exchange is clearly detected even in the present case. The qualitative tendencies of the 
phase delays are not changed by the viscous damping. 



V. CONCLUSION 

We have shown theoretically that avoided crossings can be observed in the resonance 
frequencies of acoustically coupled gas bubbles plotted as functions of the separation dis- 
tances. A state exchange taking place between the bubbles in the avoided crossing region 
has been clearly exhibited by examining the oscillation phases and transition frequencies of 
the coupled bubbles. We have clarified that the state exchange is perceived by observing the 
individual oscillation phases of the bubbles, not by observing whether the bubbles oscillate 
in-phase or out-of-phase with each other. Since the individual phase (or more properly, 
the phase difference between a bubble and the external sound) determines the sign of the 



primarily Bjerknes force |26l. l27l l28| acting on the corresponding bubble, this state exchange 
should play a role in, e.g., acoustic levitation using the force. The results of this study 
suggest that the transition frequencies introduced in Ref. [l^ can be a useful tool for de- 
tecting the state exchange, which takes place at the separation distance where an avoided 
crossing resonance frequency crosses with a transition frequency that is not a resonance 
frequency. Though we only considered triple-bubble systems in a linear arrangement, exten- 
sions to systems containing a larger number of bubbles and in different arrangements may 
be straightforward. Also, nonlinear effects on the avoided crossings and oscillation phases 
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could be examined using nonlinear models 2^, 21, 2^ 3^1 . As with other physical systems, 



the avoided crossings in acoustically coupled bubbles might be real. 
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APPENDIX A 

High-order nonlinear models for N pulsating bubbles in a liquid have been proposed, in 
which terms proportional to D^^^ (k > 2) appear that involve the translational velocities 
of the bubbles 21, 3^. In Ref. [s^l, for example, Doinikov derived the following model 



equation for spherical bubbles: 



N 



1 F- X 

-RiPi + RiPi = - — + E 1 — TvT^^i^i^j + 2RiRjR- + RiRjRAtij 
3 27ipRf \ Dl^ ^ ' ' ^ J 3 J 

R^ 

^[RiRjPj + {RiRj + bRiRj)pj\ 



2^2 ^ 



with 



R^ 

[p,- ■ (pi + 2p,-) - 3(p,- ■ ti,)[tij ■ (pi + 2p,-)]] , (A3) 



Pj - Pi 



^ f 2a\/Rio\^'' 2a AuRi ^ ,,,, 
Po + ^ ^ -^-^-Po-Pex, (A4) 

where Ri and pi are the instantaneous radius and position vector, respectively, of bubble i, Fj 
denotes external forces on bubble z, tj ^ is a unit vector, 7 is the specific heat ratio of the gas 
inside the bubbles, and /i is the viscosity. Here we showed only the incompressible version, 
though Doinikov also derived a model for bubbles in a compressible liquid. Equations 
iT| and ()A2|) represent the volume oscillation of bubble i and its translational motion. 
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respectively. The linear coupled oscillator model used in the present study is recovered from 
Eq. ()A1|) by truncating the high-order terms Hij and assuming weak driving and 7 = Kj. 

Since the velocity field forming around a pulsating sphere is proportional to where 
r is the distance measured from the center of the sphere, the truncated terms Hij, which are 
composed of the translational velocities p^, might be considered to be of fourth, or higher, 
order with respect to the inverse of the separation distances. This speculation is consistent 
with the suggestion by Harkin et al. for double-bubble systems ji^ . 

Equation ()A1|) further suggests that under the assumption of pi ~ one cannot construct 
a linear model that has higher-order accuracy than that of Eq. (^, since this assumption 
makes the high-order terms inaccurate. 



APPENDIX B 

For the convenience of readers, we show the concrete forms of |A| and CI: 
|A| 



L1L2L3 + S21S32S13 + S12S23S31 



-R10-R20-R30 

-Li(M2M3 + S23S32) - L2{MsM^ + S31S13) - L3(MiM2 + S12S21) 

+i [M1M2M3 - Mi(L2L3 - S23S32) - M2(L3Li - S31S13) 

- M3(LiL2 - S12S21)] , (Bl) 



CI = {ci,C2,C3^'^ 



[Lj - Sij){Lk - Sik) + {sij - Skj){sjk - Sik) - MjM, 



k 



+i [Mj{s,k - Lk) + Mk{sij - Lj)] 

for (z, J, k) = (1, 2, 3), (2, 3, 1), or (3, 1, 2), (B2) 



where 



For 5i ~ 0, Eqs. (jBljl and ()B2|1 reduce, respectively, to 

lAI 



L1L2L3 + S21S32S13 + S12S23S31 



-R10-R20-R30 

— L1S23S32 — L2S31S13 — L3S12S21, (B3) 
{Lj - Sij){Lk - Sik) + {sij - Skj){sjk - Sik). (B4) 



RjoRkO 

Equation ()B3|) and the real part of Eq. ()B1|) are cubic functions and Eq. ()B4|) and the real 
part of Eq. ()B2|) are quadratic functions in terms of X. The imaginary parts of Eqs. ()B1|) 
and ()B2|) can be written in a form of ujf{X), where / is quadratic in Eq. ()B1|) and linear in 
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Eq. ()B2|) . (As proved theoretically in Ref. the imaginary parts are composed of terms of 
odd orders with respect to M that are proportional to ujX^ with n being a positive integer.) 



[1] In the fields of vibration engineering and some others, avoided crossings are called "curve 
veering," "eigenvalue veering," or "(natural) frequency loci veering." They are also sometimes 
called "avoided level crossings." 

[2] J. R. Walkup, M. Dunn, D. K. Watson, and T. C. Germann, Phys. Rev. A 58, 4668 (1998). 

[3] T. Timberlake and L. E. Reichl, Phys. Rev. A 59, 2886 (1999). 

[4] S. Li and E. J. Heller, Phys. Rev. A 67, 032712 (2003). 

[5] S. Diirr, T. Volz, A. Marte, and G. Rempe, Phys. Rev. Lett. 92, 020406 (2004). 

[6] Y. Osaki, Publ. Astron. Soc. Jpn. 27, 237 (1975). 

[7] D. Gondek and J. L. Zdunik, Astron. Astrophys. 344, 117 (1999). 

[8] S. Shaik, A. loffe, A. C. Reddy and A. Pross, J. Am. Chem. Soc. 116, 262 (1994). 

[9] C. Zhu and H. Nakamura, Chem. Phys. Lett. 274, 205 (1997). 
[10] J. R. Kuttler and V. G. Sigillito, J. Sound Vib. 75, 585 (1981). 
[11] N. C. Perkins and C. D. Mote, Jr., J. Sound Vib. 106, 451 (1986). 
[12] C. Pierre, J. Sound Vib. 126, 485 (1988). 
[13] H. H. Yoo and S. H. Shin, J. Sound Vib. 212, 807 (1998). 
[14] R. S. Langley, Proc. R. Soc. Lond. A 455, 3325 (1999). 

[15] A. A. Mailybaev, O. N. Khillov, and A. P. Seyranian, J. Phys. A 38, 1723 (2005). 
[16] M. Ida, Phys. Lett. A 297, 210 (2002). 
[17] M. Ida, J. Phys. Soc. Jpn. 71, 1214 (2002). 
[18] M. Ida, Phys. Rev. E 67, 056617 (2003). 
[19] M. Ida, J. Phys. Soc. Jpn. 73, 3026 (2004). 

[20] A. Harkin, T. J. Kaper, and A. Nadim, J. Fluid Mech. 445, 377 (2001). 

[21] H. Takahha, T. Akamatsu, and S. Fujikawa, JSME Int. J. Ser. B 37, 297 (1994). 

[22] C. Feuillade, J. Acoust. Soc. Am. 98, 1178 (1995). 

[23] A. Shima, Trans. ASME, J. Basic Eng. 93, 426 (1971). 

[24] E. A. Zabolotskaya, Sov. Phys. Acoust. 30, 365 (1984). 

[25] For such a large separation distance, in practice the time delay effect due to the finite sound 
speed of the surrounding liquid should be taken into account, though we only present an 
idealized result for the sake of simplicity. 

[26] F. G. Blake, Jr., J. Acoust. Soc. Am. 21, 551 (1949). 

[27] T. J. Matula, Phil. Trans. R. Soc. Lond. A 357, 225 (1999). 

[28] I. Akhatov, R. Mettin, C. D. Ohl, U. Parlitz, and W. Lauterborn, Phys. Rev. E 55, 3747 
(1997). 

[29] R. Mettin, I. Akhatov, U. Parhtz, C. D. Ohl, and W. Lauterborn, Phys. Rev. E 56, 2924 
(1997). 

[30] A. A. Doinikov and S. T. Zavtrak, J. Acoust. Soc. Am. 99, 3849 (1996). 

[31] W. Lauterborn, T. Kurz, R. Mettin, and C. D. Ohl, Adv. Chem. Phys. 110, 295 (1999). 

[32] L. A. Crum, J. Acoust. Soc. Am. 57, 1363 (1975). 

[33] A. A. Doinikov, J. Acoust. Soc. Am. 116, 821 (2004). 



10 




FIG. 1: Resonance frequencies (U^res (rad/s) of two coupled bubbles for di ~ normalized by 
LOio (rad/s), as functions of the normalized separation distance li2- The dashed line denotes the 
monopole resonance frequency of a bubble that will be coupled with the former two bubbles in the 
next example. 
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(a) L. = 100 (b) L. = 50 (c) L. = 20 




12 ''12 ''12 



FIG. 2: Resonance frequencies Wres (rad/s) of three coupled bubbles for (5j ~ normalized by wio 
(rad/s), as functions of the normalized separation distance li2. Figures E{a), Etb), and|2fc) are for 
^23 = 100 , 50, and 20, respectively. The dashed lines denote the resonance frequencies when bubble 
3 is decoupled. 
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(a) U = 100 (b) = 50 (c) = 20 




'^12 

FIG. 4: Same as Fig. [21 but bubble 3 is smaller than the others. The bubbles are aligned as shown 
in Fig.OIb). 
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FIG. 5: Transition frequencies ojtr (rad/s) of three coupled bubbles for (5j ~ with I23 = 20 
normalized by ujiq (rad/s), as functions of the normalized separation distance li2- ujtri denotes the 
transition frequencies of bubble i. 
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FIG. 6: Phase delays 0j (rad) normalized by tt as functions of u/ujiq for different Z12 [(a), (b): for 
1x2 < hnt, (c), (d): for I12 > knt]- The solid and the dashed curves denote <pi for negligible and non- 
negligible damping, respectively, and the vertical dotted lines indicate the two lowest resonance 
frequencies (the higher is thus the second-highest resonance frequency a;2nd)- 
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